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We report on the first measurement of the differential cross section of 0-meson photoproduc- 
tion for the d("y,pK + K~)n exclusive reaction channel. The experiment was performed using a 
tagged-photon beam and the CEBAF Large Acceptance Spectrometer (CLAS) at Jefferson Lab. A 
combined analysis using data from the d(j,pK + K~)n channel and those from a previous publica- 
tion on coherent <J> production on the deuteron has been carried out to extract the <j> — N total cross 
section, a^M- The extracted 4> ~ N total cross section favors a value above 20 mb. This value is 
larger than the value extracted using vector-meson dominance models for <j> photoproduction on the 
proton. 

PACS numbers: 13.60.Le, 24.85. +p, 25.10.+S, 25.20.-x 



Multi-gluon exchange between hadrons, known as 
Pomeron exchange, is a fundamental process and plays 
an important role in high-energy interactions. At lower 
energies, this exchange manifests itself in a QCD van 
der Waals interaction |l[ . Studying multi-gluon exchange 
at lower energies is challenging because at low energies 
hadron-hadron interactions are dominated by quark ex- 
change. However, multi-gluon exchange is expected to be 
dominant in the interaction between two hadrons when 
they have no common quarks. The <p meson is unique in 
that it is nearly an ss state and hence gluon exchange is 
expected to dominate the q> — N scattering process. 

Direct measurement of the <f> — N cross section is not 
possible due to lack of <f> meson beams. An upper limit 
of (70 at ~ 11 mb 0] is obtained using the <f> photoproduc- 
tion data on the proton and the vector meson dominance 
(VMD) model [3(, which is in agreement with the esti- 
mate from the additive quark model In a geometric 
interpretation of hadron-proton total cross sections 
the radius of the <f> meson can be estimated from the 
comparison of the total cross sections of 7r — N (ovjv) and 
4> — N (ct^aO scattering. The value of ova? is ~ 24 mb [B| 
and the it radius rv is ~ 0.65 fm An upper limit of 
&<I>n ~ 11 mb leads to an upper limit value of ~ 0.43 
fm for r^. 

However, from the observed ^-dependence of nuclear <f> 
photoproduction, a larger value of (inelastic a 1 ^ as ~ 35 
mb Q, which is part of the total o^at) is obtained, which 
suggests a larger r§ value than 0.43 fm. Medium modifi- 
cation of the vector meson properties 0] (such as radius) 
or channel coupling effects [2j have been proposed to ex- 
plain the aforementioned difference in the cross section 
for the 4> meson. A similar phenomenon has been ob- 
served for the J/^i! meson 0] and the color transparency 
effect is proposed @ to explain the observation. 

In this paper, we present a determination of ct^n using 
the differential cross section of the incoherent 0-meson 
photoproduction from deuterium. This process takes ad- 
vantage of the rescattering of a 4> meson from the specta- 



tor nucleon. In the reaction 7 + d — » (fi+p + n, the rescat- 
tering process will dominate for the kinematics where 
both nucleons are energetic. The deuteron is a system 
of loosely bound nucleons and, hence, nuclear medium 
effects should not play a significant role in the <j> — N 
scattering process. In the analysis of incoherent 0-meson 
photoproduction from deuterium, the <t>—N interaction is 



parametrized as ^| oc a 2 



4>N 



111 . Il2j , where P^n 



characterizes the i-dependence of the differential cross 
section, and t is defined as the four-momentum transfer 
squared between the photon and the <f> meson. A % 2 anal- 
ysis is performed for both processes (7 + d — ► <fi + d as in 
[13j and 7 + d — > <j> + p + n) in this paper to constrain 
the values of a^N and P^n- 

The rescattering of a <f> meson off a nucleon in the 
deuteron was used to study the <p — N interaction in a 
recent analysis of CLAS glO data using coherent <f> pho- 
toproduction, 7 + d — > 4> + d In the coherent <f> 
production process, rescattering dominates in the high 
\t\ region. Results of Ref.fl3j agree with the VMD pre- 
diction (where fyff is assumed to be the one in the <fi- 
meson photoproduction from a nucleon), however, larger 
<70 N values showed better agreement with the data if one 
allowed the slope (P^n) of the ^-dependence of the cfi — N 
scattering process to be different. 

The reaction "f(d, 4>p)n was measured by detecting 
kaons from the (/>-meson decay ((f) — > K + K~ , branching 
ratio about 0.5) using the same data set as in Ref. [l3j |. 
A tagged-photon beam was generated by a 3.8-GeV elec- 
tron beam incident on a gold radiator (10 radiation 
length). The photon flux was measured by the CLAS 
photon-tagger system |l4( . The data were collected from 
a 24-cm-long liquid-deuterium target in the CLAS detec- 
tor [l5| at JLab. 

Events having the final state j + d — > K + +K~ +p + X 
were selected using a triple coincidence detection of a 
proton, a K + and, a K~ . Each particle was selected 
based on particle charge, momentum, and Time-of-Flight 
(TOF) information. The reaction d('y,pK + K~)n was 
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FIG. 1: (a) The missing mass squared distribution; (b) the 
invariant mass distribution for the K + K~ for photon energies 
between 1.65 to 3.59 GeV for the 7 + d -> K + + K~ + p + 
X process. Also shown in (b) is a fit using a Breit-Wigner 
function convoluted with the experimental resolution together 
with a flat background. 



identified in the missing mass squared distribution by 
the missing neutron shown in Fig. la. In the figure, 
the position of the neutron mass squared is shown by 
the dotted line. A ±3cr cut was employed to select the 
pK + K~(n) final state events. 

Once the reaction d(-f,pK + K~)n was identified, the 
number of 4> mesons was obtained by subtracting the 
background under the <f> peak (invariant mass spectrum 
of the K + and K~) in the ±3<r region (see Fig. lb). 
The K + K~ invariant mass distribution was fitted using 
a Breit-Wigner function convoluted with the experimen- 
tal resolution, plus a function to model the background in 
each photon energy and \t — to\ bin, where to is the mini- 
mum t value for a given photon energy. The background 
shape was assumed to be (131 ] : 



f{x) = a^Jx 2 - [2M K ) 2 + b (x 2 - (2M K ) 2 ) (x > 2M K ) 
f(x) = (x < 2M K ), (1) 

where x is the invariant mass of the K + K~, Mh is the 
kaon mass, and a and b are fitting parameters. The back- 
ground was also fitted with a linear shape. The results 
from fitting these two shapes were compared to estimate 
the systematic uncertainties due to the subtraction of the 
background. 

A Monte Carlo (MC) 16] simulation of the CLAS de- 
tector was carried out to determine the efficiency for the 
detection of the 7 + d — ► p + K + + K~ reaction. Two 
event generators were used in two different missing neu- 
tron momentum regions. A quasi-free event generator for 
the 7 + p — > <j) + p process in the deuteron was employed 
for the missing neutron momentum distribution below 
0.18 GeV/c, where the agreement between the MC simu- 
lation and the data is very good. The deuteron wave 
function was based on the Bonn potential 
differential cross sections from CLAS UM for 



neutron momentum greater than 0.18 GeV/c, the gener- 
ated events were weighted by dp d d £ dn from Ref.fioj 
to include both the <f> — N and N — N final state inter- 
actions. In both cases, the <f> mass distribution was mod- 
eled by a Breit-Wigner shape with the resonance centered 
at the mass of 1.019 GcV/c 2 and with a FWHM of 
0.0044 GeV/c 2 . The <^-meson decay angular distribution 
was taken as UM: 



W(cos(6 H )) 



3,1 
2 V 2 



17J. The 

photo- 
production from the proton were used. For the missing 



( 7; ( 1 -Poo ) sin2 °h +P00 cos 2 H ) +a cos 6 H , 

(2) 

where Pq is the spin density matrix element, Oh denotes 
the polar angle of the K + in the 0-meson rest frame, and 
a accounts for an interference between the <f> and the non- 
resonant S*-wave K + K~ pair production [2(|. Helicity 
conserving amplitudes give pg = 0, while single-helicity 
flip amplitudes require pg 7^ 0. A value of 0.1 (0.05) was 
used for pg (a) in the MC simulation. The MC gen- 
erated events were used as input to the GEANT3-bascd 
CLAS simulation [lj| . They were then reconstructed us- 
ing the same event reconstruction algorithm as was used 
for the data. The acceptance was obtained by the ra- 
tio of the number of events that passed the analysis cuts 
to the number of generated 4> events. The average dif- 
ferential cross section for each photon energy and \t — £q| 
bin was extracted by dividing the normalized yield (num- 
ber of selected events divided by the integrated photon 
flux including the DAQ dead time, the target thickness, 
the <p decay branching ratio and \t — t \ bin size) by the 
acceptance which includes the detection efficiency. The 
differential cross sections were then bin-centered at fixed 
t values and a finite binning correction was applied. 

Several sources contribute to the overall systematic un- 
certainty in extracting the differential cross section. The 
systematic uncertainties associated with particle identi- 
fication and the missing mass cut were determined to be 
0.5 - 7.0% and 0.5 - 5.0%, which are the values found 
across the different bins of photon energy and t, re- 
spectively These were determined by varying the cor- 
responding cuts by ±10%. The uncertainties in the pa- 
rameters of the 4> decay angular distribution, Pq Q and a, 
were 10% and 5% [3, respectively, leading to 1.5- 
6.5% systematic uncertainties. The background obtained 
from the non-linear background shape was on average 
5% smaller than that from the linear background. The 
systematic uncertainties from the acceptance dependence 
on the cross section model varied from 0.5% to 9%. The 
uncertainty in the photon flux was 5% 21 , 22| ■ The un- 
certainty of the bin-centering corrections were typically 
between 0.5% and 6.0%, based on our current knowl- 
edge of the CLAS detector and the uncertainty in the 
input cross section. Combining all systematic uncertain- 
ties in quadrature, the overall systematic uncertainties 
vary from 7% — 18% depending on the kinematics. 

In Fig. [21 differential cross sections, ^ from this work 
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(red solid circles) for the reaction 7 + d — > <p + p + n 
are presented. In Fig. 2a and Fig. 2b, we present 4? as 
a function of |f — t | f° r a photon energy range of 1.65 
- 2.62 GeV and 2.62 - 3.59 GcV, respectively (same as 
those in Ref . [l3| ) . The detected proton and the missing 
neutron momentum span a range of 0.18 to 2.0 GeV/c 
(Figs. 2a and 2b). For the low missing momentum re- 
gion (the momenta of the reconstructed neutron smaller 
than 0.18 GcV/c), cross sections over the same photon 
energy ranges are shown in Fig. 2c and Fig. 2d. The er- 
ror bars shown are the statistical uncertainties, and the 
overall systematic uncertainties are shown by the black 
band. Also plotted are the predictions [lj| for quasi-free 
cj> production and rescattering for four sets of ct^tv an d 
fltfrN- The calculations are performed based on the model 
assumptions of pomeron exchange or the two-gluon ex- 
change interaction in the cf> — N rescattering process. 

In [10( , the neutron-proton rescattering amplitude has 
been taken into account and treated in the same way as in 
the analysis of the d(e, e'p)n channel 23l.l24|. The models 
for the 6 photoproduction on the proton [25l. I2H ] used in 
Ref. II 011 describe the published experimental cross sec- 
tions [20|, [27 1 reasonably well for photon energies in the 
vicinity of 3.4 GeV and above. However, at low photon 



energies the model underestimates the experimental cross 
sections of 7+p — ► 4>+p [l8| above \t — to| = 1 (GeV/c) 2 
by a factor that can reach 10 at 2.5 (GeV/c) 2 . The ex- 
tracted cross sections [28| of the 7 + "p" — > p + (f> from 
the d(j, 4>p)n channel at low neutron missing momentum 
from this work are in good agreement with the extracted 
cross sections of 7 + p ^ cf> + p pj| . The flattening be- 
havior of the experimental cross sections with increasing 
\t — to\ values is well accounted for by the coupling be- 
tween the 4> and uj production channels [2^, which is 
discussed later. This channel coupling effect has not yet 
been implemented in Laget's code that we use, and its 
effect will be investigated in a future study. 

We simply note that in the high missing momentum re- 
gion where contributions of rescattering processes {(p — N 
or N — N) are significant, the dominant contribution to 
the matrix element comes from the photoproduction of <f> 
on a nucleon at rest, which gets its measured momentum 
in the rescattering process (wc refer to Ref. [l(J for a de- 
tailed discussion). The consequence is that in the scat- 
tering loop the elementary photoproduction amplitude 
~f + N^<fr + N is almost the same as in the quasifree 
case, when the spectator is at rest. The amplitude for 
the j + d—^c^+p + n process can be written as: 



J 



A = A^ ® (1 + + Agf) + Ajf ® (1 + Ag + A%) = {Af p + A«f) 0(1 + Ag + A%) 



(3) 



where A® N , A^ and A^ are the amplitude for the (/>- 
meson photoproduction from nucleon, the amplitude for 
the proton-neutron final state interaction, and the am- 
plitude for the 4>-N final state interaction, respectively. 
The second step assumes isospin symmetry for the (fi-N 
interaction. The <X> represents a convoluted integral over 
internal momentum. Thus, deviations from data in the 
model cross section of <j) photoproduction from nucleons 
will lead to deviations of the calculated cross sections in 
the high spectator nucleon momentum region from the 
data, even without final state interaction effects. In or- 
der to minimize this effect, we form the following ratio 
between the results with a high spectator momentum cut 
and the results with a low spectator momentum cut for 
each value of It — to\: 



f 2GeV _da 
j£ _ °H_ _ JO.lSGeV dtdPrr. 



■dP„ 



f 

Jo 



0.18GeV da 



GeV dtdP, r 



■dP„ 



(4) 



The amplitude with a low spectator nucleon momentum 
cut is 



a _ aQS 1 aQS 



(5) 



within the quasi-free reaction mechanism. Therefore, the 



extracted ratio is 



R 



Uf p +Af n )®{\ 



A fsi 

pn 



U Q " + A^)\ 



(G) 



i ( i + A£r+A^)i 5 



The second step relies on the factorization approxima- 
tion which works better when the elementary amplitude 
varys slowly with energy. This is the case for the 0-meson 
photoproduction channel. This implies that to first order 
in the ratio R, the elementary (f> photoproduction ampli- 
tude and, hence, the related model uncertainties, should 
cancel out. 

The comparison of the experimental ratios to the 
model calculations from Ref. [lfj are presented in Fig. 
3. The experimental data are shown with statistical un- 
certainties. The systematic uncertainties are shown by 
the black band. In the low energy range, E 1 = 1.65-2.62 
GeV (Fig. 3a), the data are best described by the pa- 
rameters of (70AT = 10 mb and ft^N — 6 (GeV/c)~ 2 . In 
the range E 1 = 2.62-3.59 GeV, the data can be described 
well by all four calculations shown in Fig. 3 including the 
rescattering effect. 

In order to constrain the value of o^n using our data, 
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FIG. 2: Differential cross section ^ from ~/(d, K + K~ p)n pro- 
cess for photon energies of 1.65-2.62 GeV (a, c) and 2.62- 
3.59 GeV (b, d). The missing momentum is higher than 0.18 
GeV/c in (a) and (b), and lower than 0.18 GeV/c in (c) and 
(d). The results of this work are shown in red solid circles. 
The black bands represent the systematic uncertainties. The 
label "30 + 10" indicates the calculation from Laget [l(j with 
vtot = 30 mb and /? 0JV = 10 (GeV/c)" 2 . The legend for the 
calculations is presented for better visibility and it applies to 
all four panels. 



a x 2 analysis was performed for results from both the 
-f + d^<fi+p + (n) process (R = and the published 
7+d — ► 4>+d coherent channel (l3j by mapping the phase 
space of a^N and fl^N- The x 2 is defined as: 



da t a — Real) 



i=l 



5R data 



(<) 



where N is the number of data points. Rdata (Real) is the 
cross section ratio defined in Eq. 3 for data (calculation) . 
The SRdata is quadrature sum of the point-to-point sys- 
tematic uncertainty and the statistical uncertainty. For 
the coherent process, Rdata (Real) is the differential cross 
section data (calculation) . The calculations of Laget [l(| , 
which include pomcron exchange for the elementary cj> 
meson photoproduction cross section on the nucleon, are 
used for the ^y + d—tifi + p + n channel. The pomeron 
and the two-gluon exchange versions of the model lead 
to very similar results at \t — to I values smaller than 1.5 
(GeV/c) 2 , as can be seen in Fig. 2 and Fig. 3. For the 
coherent pro duction channel, calculations from Sargsian 
et al. [ill. Il2l] are used. 
Fig. 4a and Fig. 4b show the confidence level for both 



processes in the two photon energy regions. While the 
energy dependence in g^n and f3<pN might be a possible 
explanation for the difference between Fig. 4a and 4b, 
the combined analysis favors a value of o^at larger than 
20 mb. Our results are consistent with that extracted 
from the SPring-8 data Q in which Li, C, Al and Cu 
nuclear targets were used. Further more, our combined 
analysis gives a lower bound of 6 (GeV/c) -2 for the j3<f,N 
parameter. 

Medium modifications that have been suggested by the 
SPring-8 data @ can hardly explain a large g^n cross 
section (large r^) in deuterium. More likely it may re- 
flect the fact that other mechanisms, beyond np and 4>N 
rescattering, are at play and are more important than 
the medium modifications, for example the QCD van der 
Waals force. On the one hand, the coupling of the 4> (via 
two-gluon exchange) to hidden color components [30j in- 
side the deuteron may contribute to large missing mo- 
menta and leave less room for <p — N rescattering. On 
the other hand, the coupling to a cryptoexotic baryon 
(baryon with hidden strangeness) , = uddss, may also 
contribute Q. However, the most likely explanation lies 
in the to — <p mixing. The photon produces an u) meson 
on one nucleon, which elastically scatters on the second 
nucleon before transforming into a 4> meson. The corre- 
sponding matrix element has the same structure as the 
clastic <fi rescattering matrix element that we considered 
in our analysis. An effective u^n cross section value can 
be written as: 



eff 
'4>N 



l(T 7 N->uN 











(8) 



50 [13, [3l| in 



<7 7 jv_.0jv 

a —t region of 1 to 3 (GeV/c) 2 , and the LO-<p mixing co- 
efficient g^^^ ~ 0.09 [32j, one can reconcile the effective 
<pN cross section of ~ 30 mb with the VMD values of the 
<Ju>n (~ 25 mb) and ct^n (~ 10 mb) [2{|. One way to put 
the (j> — N cross section on more solid ground would be to 
select the part of the phase space where the <f> — N rescat- 
tering dominates, using the method proposed in Ref. [Io| . 
Future high statistics data from a luminosity upgraded 
CLAS12 [33j detector will help disentangle these possibil- 
ities. Furthermore, future improved and new theoretical 
calculations will allow us to study the model uncertainty 
in the extraction of the <^>-N total cross section. 
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FIG. 3: Cross section ratio between the high and the low missing momentum regions for photon energies of 1.65-2.62 GeV 
(a) and 2.62-3.59 GeV (b). The results of this work are shown in red solid circles. The black bands represent the systematic 
uncertainties. We use same notations as those in Fig. [2] for calculations from Laget 
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FIG. 4: The 70% (shaded area) and the 95% (open area) confidence level plots shown for the 7 + cj> + p + (n) channel 
(red), the 7 + d — + (f> + d coherent channel (black) for photon energies of 1.65-2.62 GeV (a) and 2.62-3.59 GeV (b). 
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